Structure of Cubic Matrix Mechanics 

Yoshiharu Kawamura 1 
Department of Physics, Shinshu University, Matsumoto 390-8621, Japan 

Abstract 

We study the structure of cubic matrix mechanics based on three- 
index objects. It is shown that there exists a counterpart of canonical 
structure in classical mechanics. 
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1 Introduction 



The study of the basic structure of mechanics is especially significant for 
several reasons. First, such a study provides a unified description applicable 
to various systems. For example, the canonical formalism in classical me- 
chanics (CM) has been applied to a wide range of dynamical systems, from 
a system of point particles to continuous systems, such as relativistic field 
theories. [1] Second it is useful in making clear the meaning of physical quan- 
tities and symmetries of the system. Third it is informative for the purpose 
of constructing a new mechanics based on the structure of established me- 
chanics. In fact, quantum mechanics (QM) was constructed with the aid of 
experimental results and the correspondence principle. [2] 

If there is a new mechanics beyond QM, it is natural to suppose that it 
must possess a structure similar to that of QM and CM, because a new scheme 
would contain QM as a limiting case, and the correspondence principle relates 
QM to CM. Based on this conjecture, we consider the following questions: 

1. What is the basic structure of such a new mechanics? 

2. Is there a new, generalized mechanics beyond QM and CM? 

The first question can be rephrased as a question of what structure of 
QM and CM is preserved and how other parts are modified. We expect 
that this new scheme would contain descriptions of time development and 
symmetry transformation properties of a system that are similar to those 
of QM and CM. therefore, we require that the algebraic structure of equa- 
tions of motion and symmetry transformations be preserved ( up to anomalous 
breakings) in the new mechanics. Regarding the second question, a general- 
ization of Heisenberg's matrix mechanics based on many-index objects has 
been proposed. [3] It is important to explore the entire structure in order to 
obtain information on the physical meaning and relation to reality. 

In this paper, we study the basic structure of cubic matrix mechanics and 
show that this mechanics possesses a counterpart to the canonical structure 
in CM. 

This paper is organized as follows. In the next section, we review the 
canonical structure of classical mechanics and discuss the basic structure 
that a mechanics beyond CM should possess. The structure of cubic matrix 
mechanics is studied in 53. Section 4 is devoted to conclusions and discussion. 
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2 Classical mechanics and beyond 



2.1 Canonical structure of classical mechanics 

Here we review the canonical structure of classical mechanics.fi] The canoni- 
cal variables q a = q a (t) and p a = p a (t) satisfy Hamilton's canonical equations 

dq a _ dH_ dpa _ 
dt dpj dt dq^ 

where a = 1, • • • ,n and H is the Hamiltonian. Physical variables are given 
by functions of the canonical variables and the time variable t; e.g., A = 
A(q a ,p a ,t) and B = B(q a ,p a ,t). Hereafter we consider systems such that 
physical variables do not contain t explicitly, that is, closed physical systems. 
The Poisson bracket of two variables A and B with respect to q a and p a is 
defined by 

» [OA 8B OA dB\ , o . 
Hence, the Poisson brackets of the canonical variables are given by 

{q a ,qb}pB = {Pa,Pb}pB = 0, {q a ,Pb}pB = Sab. (3) 

The basic features of the Poisson bracket are as follows: 

{A,B} PB = -{B,A} PB , (antisymmetry) (4) 

{A + B, C} PB = {A, C} PB + {B, C} PB , (linearity) (5) 

{{A, B} PB , C} PB + {{B, C} PB , A} PB + {{C, A} PB , B} PB = 0, (6) 

{AB, C} PB = A{B, C} PB + {A, C} PB B. (derivation rule) (7) 

By use of (1) and (2), the physical variable A is shown to satisfy the equation 

d 4 = {A,HW (8, 

A transformation A — > A' = C(A) that preserves the Poisson bracket 
structure is called 'canonical': 

{A, B} PB — C({A, B} PB ) = {C(A),C(B)} PB . (9) 
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The infinitesimal version A — > A' = A + 5 A is given by 

SA = {A,G} PB Ss, (10) 

where G is the generator of the transformation and 5s is an infinitesimal 
parameter. We can show that the Poisson bracket structure is preserved un- 
der the transformation (10) by using the Jacobi identity (6). When C(H) = 
H, the transformation is a symmetry transformation, and its generator G 
is a constant of motion, i.e., dG/dt = {G,H} PB = 0. Then the equa- 
tion of motion is form-invariant under the canonical transformation, so that 
dC(A)/dt = {C(A),C(H)} PB . If Gj are conserved quantities, {Gi,Gj}p B are 
also conserved, as seen from (6) and (8). 

2.2 Beyond classical mechanics 

The structure of classical mechanics is so simple and elegant that we expect 
that a mechanics beyond CM must have a similar structure. In this subsec- 
tion, we present a conjecture for the basic structure of a new mechanics by 
analogy to the canonical structure of CM.t 

First we require the new scheme to possess the following properties: 

1. There are counterparts of the canonical variables in CM, which are de- 
noted Q a = Q a (t) and P a = P a (t), and physical variables are functions 
of Q a and P a in a closed system. * There exists a counterpart of the 
Poisson bracket, which we call 'generalized bracket', and the bracket 
relations for Q a and P a are conditions that place restrictions on the 
phase space. The generalized bracket does not necessarily possess all 
the algebraic properties of the Poisson bracket. However, at least they 
possess properties of antisymmetry and linearity. 

2. The same type of equation of motion holds for physical variables. More 
specifically, the equation of motion is obtained by the replacement of 
the Poisson bracket with the generalized bracket. 

^Though QM has been established as a realistic mechanics beyond CM, we temporarily 
set QM aside to give full freedom to our imagination in the construction of a new scheme 
in this subsection. 

■^In this paper, we do not consider a system with extented phase space variables, de- 
scribed by Nambu.[4] It is an interesting subject to study the relation between Nambu 
mechanics and generalized matrices. [5] 
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3. There is a transformation that preserves the generalized bracket struc- 
ture that we call 'a generalized canonical transformation'. The Jacobi 
identity for the generalized bracket holds if it contains a conserved 
quantity. Continuous symmetry transformations are realized as gener- 
alized canonical transformations generated by conserved quantities. 

Next, we formulate the basic structure of a new mechanics based on the 
above requirements. 

1. Let us denote the generalized bracket by £>(*,*) and impose the fol- 
lowing conditions on Q a and P a : 

B(Q a , Q b ) = B(P a , P b ) = 0, B(Q a , P b ) = 5 ab Q, (11) 

where is a constant of motion, and the bracket of and a conserved 
quantity A vanishes, i.e. £>(0, A) = 0. The antisymmetry and linearity 
conditions are expressed by 

B(A,B) = -B(B,A), (12) 
B(A + B,C) =B(A,C) +B(B,C). (13) 

We do not necessarily require the Jacobi identity nor the derivation 
rule as a property of £>(*, *) for generic variables. 

2. The equation of motion for a physical quantity A is given by 

HA 

1S = B(A,H), (14) 

where H corresponds to the Hamiltonian and is interpreted as the gen- 
erator of time evolution from the symmetry property discussed just 
below. 

3. A generalized canonical transformation is defined by the transformation 
A — > A' = G(A), which preserves the structure of £>(*, *): 

B(A, B) — g(B(A, B)) = B(Q(A), 0(B)). (15) 

The infinitesimal version of (15) is written 

SB(A,B)=B{SA,B)+B(A,SB), (16) 
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under the infinitesimal generalized canonical transformation A — > A' = 
A + 8 A. For a conserved quantity G, i.e., dG/dt = B(G,H) = 0, the 
following Jacobi identity holds: 

B(B(A, B),G) + B(B(B, G),A) + B(B(G, A),B) = 0. (17) 

Then, a symmetry transformation is given by the infinitesimal gener- 
alized canonical transformation,^ 

8 A = B(A, G)6s. (18) 

For conserved quantities Gi, B(Gi, Gj) are also conserved quantities, as 
seen from (14) and (17). 

2.3 Canonical formalism of quantum mechanics 

It is well known that there is a correspondence between the algebraic struc- 
ture of quantum mechanics and that of classical mechanics. [6] Here we give a 
brief review of the canonical formalism of quantum mechanics using Heisen- 
berg's matrix description for later convenience. Physical quantities are rep- 
resented by hermitian matrices that can be written 

A mn (t) = A mn ^ = A mn e^- E ^\ (19) 

where the phase factor implies that a change in energy E m — E n appears 
as radiation with angular frequency Q mn , and the hermiticity of A mn (t) is 
expressed by A* nm {t) = A mn (t). The quantity A mn (t) is a function of canon- 
ical pairs {Q a )mnif) and (P a ) m „(t). A generalized bracket is a commutator 
divided by ih: 

B(A, B) mn = lj2( A m k(t)B kn (t) - B mk (t)A kn (t)) = ^[A(t),B(t)] mn . (20) 

k 

The commutators for the canonical pairs are given by 

[Qa(t), Qb{t)}mn = [P a (t) , P b (t)] mn = 0, [Q a (t), P b (t)] mn = lh8 ab 8 mn . (21) 

§It is not certain whether every continuous generalized canonical transformation A — > 
A' = Q{A) can be constructed from the infinitesimal one given by 8 A = B(A, F)8s, where 
F is a generator. Here, we require the algebraic structure of symmetry transformations to 
be identical to that in CM. 
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It is obvious that the generalized bracket (20) possesses every property, in- 
cluding the Jacobi identity and the derivation rule, possessed by the Poisson 
bracket. 

A physical quantity A mn (t) satisfies the Heisenberg equation 

^-A mn (t) = l -{E m - E n )A mn {t) = ±r[A(t),H] mn , (22) 
at h in 

where the Hamiltonian H is a diagonal matrix written H mn = E m 5 mn . 

A generalized canonical transformation is a unitary transformation given 

by 

A mn (t) - A' mn (t) = (rfAU) mn (t), (23) 
where U mn is a unitary matrix. The infinitesimal version is given by 

5A mn (t) = ^-[A(t),G] mn 5s. (24) 

For a unitary transformation generated by conserved quantities (Gi) mn , [Gi, Gj] 
are also conserved. 



3 Cubic matrix mechanics 

We have discussed the basic structure that a new mechanics beyond classical 
mechanics should possess. Quantum mechanics is a typical example, and 
it describes the microscopic world very successfully, but there is no definite 
reason that QM is the unique mechanics capable of describing nature at 
a fundamental level (around and beyond the gravitational scale). For this 
reason, it is still meaningful to construct a new, generalized mechanics and 
study its properties. In this section, we study the structure of cubic matrix 
mechanics, which has been recently proposed. [3] 

3.1 Cubic matrix 

Here we state our definition of a cubic matrix and its related terminology. A 
cubic matrix is an object with three indices, Ai mn , which is a generalization 
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of a usual matrix, such as B mn .* We refer to a cubic matrix whose elements 
possess cyclic symmetry, i.e., A\ mn = A mn \ = A n i m , as a cyclic cubic matrix. 
We define the hermiticity of a cubic matrix by A Vm i n i(t) = A* lmn (t) for odd 
permutations among indices and refer to a cubic matrix with hermicity as 
a hermitian cubic matrix. A hermitian cubic matrix is a special type of 
cyclic cubic matrix, because there it obeys a relation A[ mn = A* mln = A mn i = 
A* nml = A nim = A* lnrn . We refer to the following form of a cubic matrix as a 
normal form or a normal cubic matrix: 

d n l + dnldlm- (25) 

A normal cubic matrix is also a special type of cyclic cubic matrix. The 
elements of a cubic matrix are treated as c-numbers throughout this paper. 

3.2 Cubic matrix mechanics and its structure 

The basic ingredient of this mechanics is a cyclic cubic matrix given by 

Aimnit) = A\ mn e lmn , (26) 
where the angular frequency fli mn has the form 

Qlmn = Wlm ~ ^In + ^mn = (5w)lnm, ^ml = —Wlm- (27) 

We assume the generalization of Bohrs' frequency condition" 

Qlmn = J^Elm ~ E in + E mn ) , (28) 

where E[ m (= —E m i) are energy eigenvalues. The angular frequencies fl[ mn 
have the properties 

^I'm'n' = Sgn(P)f2 imn , (SQ)[ mn k = Qlmn — ^Imk + ^Ink ~ ^mnk = 0, (29) 

^Awata, Li, Minic and Yoneya introduced many-index objects to quantize Nambu 
brackets. [7] We find that our definition of the triple product among cubic matrices is 
different from theirs, because we require a generalization of the Ritz rule in the phase 
factor, but not necessarily the associativity of the products. 

II Here and hereafter we use the reduced Planck constant h = as the unit of action 
in our cubic matrix mechanics for simplicity. 
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where sgn(P) is +1 and —1 for even and odd permutation among indices, re- 
spectively. The operator 5 is regarded as a coboundary operator that changes 
k-th antisymmetric objects into (/c+l)-th objects, and this operation is nilpo- 
tent, i.e. S 2 (*) = 0.[8] The frequency f^ mn is regarded as a 3-coboundary. 

If we define the triple product among cubic matrices Ai mn (t) = Ai mn e lillmnt , 
B lmn (t) = B lmn ^ n t and Clmn (^ = C lmn e in ^nt by 

(A(t)B(t)C(t)) lmn = ]T A lmk (t)B lkn (t)C kmn (t) = {ABC) lmn ^™\ (30) 

k 

this product takes the same form as (26) with the relation (29), which is 
a generalization of the Ritz rule. We comment that the resultant three- 
index object (ABC)i mn e %nimnt does not always have cyclic symmetry, even if 
Aimnit)) Bimnit) and Ci mn (t) are cyclic cubic matrices. Note that this product 
is, in general, neither commutative nor associative; that is, (ABC)i mn ^ 
(BAC) lmn and (AB(CDE)) lmn ^ (A(BCD)E) lmn ^ {{ABC)DE) lmn . The 
triple-commutator is defined by 

[A(t), B(t), C(t)] lmn = (A(t)B(t)C(t) + B(t)C(t)A(t) + C(t)A(t)B(t) 
- B(t)A(t)C(t) - A(t)C(t)B(t) - C(t)B(t)A(t)) lmn . (31) 

If Ai mn (t), B lmn (t) and Ci mn (t) are hermitian matrices, i[A(t), B(t), C(t)]i mn 
is also a hermitian cubic matrix. The generalized bracket is defined by use 
of the triple-commutator (31) as 

B(A,B) lmn ee hA(t),I,B(t)] lmn , (32) 
in 

where / is a special type of normal cubic matrix given by 

hmn = ^m(l — $mn) + ^mn(l — $nl) + $rd(l ~ $lm)- (33) 

By definition, we find that the generalized bracket (32) has the properties of 
antisymmetry and linearity, as seen from the relations 

[A(t), I, B(t)] lmn = ~[B(t), I, A(t)] lmn , (34) 
[A(t) + B(t), I, C(t)] lmn = [A(t), I, C(t)] lmn + [B(t), /, C{t)] lmn . (35) 

Note that neither the derivation rule nor the Jacobi identity necessarily holds 
for generic variables. (See the appendix for features of the triple-commutator 
[A, I, B}.) 
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We impose the following conditions on the canonical pairs (Q a )i mn and 

(yPajlmn- 

[Q a (t),I,Q b (t)] lmn = [P a (t),I,P b (t)] lmn = 0, 

[Qa(t), I, Pb(t)]lmn = ^ab^lmn- (36) 

Here Qi mn can be a normal cubic matrix, because the conditions should be 
time-independent, and a normal cubic matrix is a constant of motion, as seen 
below. 

The cyclic cubic matrix A[ mn (t) yields the generalization of the Heisenberg 
equation, 

-^Aimn(t) = iVti mn Ai mn (t) = —[A(t),I,H]i mn , (37) 
where H is the Hamiltonian written 

Hlmn = -blmEmn + -&mnEnl + -&nlEl m - (38) 

Because the Hamiltonian is a normal form, we find that an arbitrary normal 
cubic matrix is a constant of motion: ih(dA^ /dt)i mn = [A^ N \ I, H]i mn = 
0. 

The generalized bracket structure (32) is preserved by the infinitesimal 
transformation 

SA lmn (t) = ^[A(t),I,G^] lmn Ss, (39) 

where is a normal cubic matrix, that is, a generator of the symmetry 
transformation. Here, we use the fact that the Jacobi identity holds for G^ N \ 
so that 

[[A, I, B], I, G^] lmn + [[B, I, GW], /, A] lmn + [[(?("> /, A], I, B] lmn = 0.(40) 

We find that the derivation rule 

[ABC, I, G^U = (AB[C, I, G^]) lmn + (A[B, I, G^]C) lmn 

+ ([A,I,G (N ^]BC) lmn (41) 

holds for G^ if (ABC)i mn is a cyclic cubic matrix. 
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4 Conclusions and discussion 



We have studied the basic structure of a cubic matrix mechanics and shown 
that this mechanics meets the requirement that the algebraic structure of 
equations of motion and symmetry transformations be preserved (up to anoma- 
lous breakings) in any new mechanics beyond quantum mechanics and clas- 
sical mechanics. 

The basic structure of this mechanics is summarized as follows. There is 
a symmetry transformation for a physical quantity A [which is a cyclic cubic 
matrix in the cubic matrix mechanics] , whose infinitesimal version is given by 
5 A = B(A, G)5s, where £>(*, *) is the counterpart of the Poisson bracket in 
CM [B{A, G) = jp-[A, I, G] in cubic matrix mechanics], and G is the generator 
of the transformation. The time evolution of A is regarded as a symmetry 
transformation generated by the Hamilotonian H: 5 A = B(A, H)5t [which is 
the generalization of the Heisenberg equation (37) in cubic matrix mechanics] . 
The generator G is a constant of motion, i.e. dG/dt = B(G,H) = 0. [The 
normal form G^ n is a conserved quantity in cubic matrix mechanics]. The 
Jacobi identity holds in the case that it contains a conserved quantity such 
as (17) [(40) in cubic matrix mechanics]. The bracket structure is preserved 
under the symmetry transformation generated by a normal cubic matrix in 
cubic matrix mechanics, as seen from the Jacobi identity. 

The above-stated requirement stems from the expectation that a conser- 
vation law survives beyond QM, but it may be too naive and need some 
modifications if we wish to incorporate a gravitational interaction. In fact, 
there is the conjecture that no continuous global symmetry exists in a quan- 
tum theory including gravity. [9] Moreover, the theory should be formulated 
in a background-independent way, as the theory of general relativity. There- 
fore, the new scheme discussed in this paper can be interpreted as an effective 
description of an underlying mechanics after fixing the background geometry 
and ignoring dynamical degrees of freedom for the graviton. 

There are interesting subjects that remain. One is the construction of a 
matrix mechanics in which physical variables are n-index objects (n > 4). 
This was proposed in Ref . [3] , but its formulation has not yet been completed. 
A conjecture for "Hamiltonians" and a generalized bracket has been given, 
but there is the limitation that the generalized Heisenberg equation holds 
only for n-index objects whose indices are completely different. Another 
subject is the study of the relation between cubic matrix mechanics and 
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Heisenberg matrix mechanics. Such studies are now in progress. 
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In this appendix, we study features of the triple-commutator [A, I, B] for 
cyclic cubic matrices Ai mn and Bi mn . This commutator is written 

[A, I, B]i mn = Ai mn B[ mn — Bi mn Ai mn + (AB)® mn , (42) 

where Bi mn and (AB)® mn are defined by 

Bimn = Bi mm B m n -\- B mnn B nmm -\- B n n B\ nn (43) 

and 

(AB)® mn = 5[ m '^ / (A mn i,B nrn i i — B m nkA n mk) + 5 mn ^(^4 n ifc-B/nfc ~~ B n i k A[ n k) 

k k 

+ 8 n i y^{A lmk B m i k — Bi mk A m i k ), (44) 
k 

respectively. The features of Bi mn and (AB)° mn are as follows: 

1. Bimn is skew-symmetric: 

Bimn B mn i B n [ m B nm [ B m i n B[ nm . (45) 

2. Bi rnn is a 5-closed (3-cocycle): 

(SB)i mnk = Bimn — Bi nl k + Anfe — B mn k = 0. (46) 

3- {AB)1 mn is a normal cubic matrix with (AB)? mn = -(BA)? mn . 
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For a normal cubic matrix A\^ n = 5i m a mn + 5 mn a nl + 5 nl a lm , A( N ) lmn is 
given by 

A&hmn = A^\ mm -A^ mU + A^ mnn -A^ nmm + A^ nll -A^ lnn 
= ^(ati + a^ + atf), (47) 

where a{~ } = |(a Jm -a ml ). 

We obtain the following relations for cyclic cubic matrices from the above 
expressions and features. 

1. For an arbitrary cyclic cubic matrix Ai mn and an arbitrary normal 
cubic matrix the relations {AB^f lmn = and [A, I, B^ N \ mn = 
A lmn BW lmn hold. 

2. Arbitrary normal cubic matrices A\^ n and B^l commute; 

[^,/,5 ( \ n = 0. (48) 

3. The Jacobi identity holds if any of A, B or C has a normal form seen 
from the relation 

[[A, I, B], I Mi™ + [[B, I, ClI, A} lmn + [[CJ^A], J, B} lmn 
= -(AB) lrnn C lmn - (BC)° lmn A lmn - (CA)° lmn B lmn .(49) 

4. The derivation rule holds for arbitrary normal cubic matrices -D^2, so 
that 

[ABC, I, DW] lmn = (AB[C, I, DW]) lmn + (A[B, I, D^]C) lmn 
+ ([A,I,D (N ^]BC) lmn , (50) 

if (ABC)i mn is a cyclic cubic matrix. Here we have used the fact that 
D( N \ mn is a 3-cocycle. 
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